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Abstract 

We introduce the general arbitrage-free valuation framework for counterparty risk 
adjustments in presence of bilateral default risk, including default of the investor. We 
illustrate the symmetry in the valuation and show that the adjustment involves a long 
position in a put option plus a short position in a call option, both with zero strike 
and written on the residual net value of the contract at the relevant default times. 
We allow for correlation between the default times of the investor, counterparty and 
underlying portfolio risk factors. We use arbitrage-free stochastic dynamical models. 
We then specialize our analysis to Credit Default Swaps (CDS) as underlying portfolio, 
generalizing the work of Brigo and Chourdakis (2008) |10j who deal with unilateral 
and asymmetric counterparty risk. We introduce stochastic intensity models and a 
trivariate copula function on the default times exponential variables to model default 
dependence. Similarly to [10], we find that both default correlation and credit spread 
volatilities have a relevant and structured impact on the adjustment. Differently from 
[To] , the two parties will now agree on the credit valuation adjustment. We study a 
case involving British Airways, Lehman Brothers and Royal Dutch Shell, illustrating 
the bilateral adjustments in concrete crisis situations. 
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1 Introduction 
General motivation. 

In general, the reason to introduce counterparty risk when evaluating a contract is linked to 
the fact that many financial contracts are traded over the counter, so that the credit quality 
of the counterparty and of the investor may be important. This has become even more 
relevant in recent years as some protection sellers such as mono-line insurers or investment 
banks have witnessed increasing default probabilities or even default events, the case of 
Lehman Brothers being a clear example. 

When investing in default risky assets, one requires a risk premium as a reward for 
assuming the default risk. If one thinks, for example, of a corporate bond, it is known that 
the yield is higher than the corresponding yield of an equivalent treasury bond, and this 
difference is usually called credit spread [12]. The (positive) credit spread implies a lower 
price for the bond when compared to default free bonds. Many works have been proposed 
recently to explain the term structure of credit spreads, such as [7], [S], [21], [E] and [2J, who 
also focuses on hedging. This reduction in value is a typical feature: the value of a generic 
claim traded with a counterparty subject to default risk is always smaller than the value of 
the same claim traded with a counterparty having a null default probability, as was shown 
formally for example in 

Bilateral risk and symmetry. 

This paper introduces a general arbitrage-free valuation framework for bilateral counterparty 
default risk. By 'bilateral' we intend to point out that the default of the investor is included 
into the framework, contrary to earlier works. This brings about symmetry, so that the 
price of the position including counterparty risk to the investor is exactly the opposite of 
the price of the position to the counterparty. This is clearly not the case if each of the two 
parties computes the present value assuming itself to be default-free and allowing for default 
of the other party only. This asymmetry would not matter in situations where financial 
investors had high credit quality and counterparties rather low one. Indeed, in such a case 
both parties would consider the investor as default-free and the counterparty as defaultable, 
so that inclusion of the investor default would be pointless, given that it happens in almost 
no scenario. However, recent events show that it is no longer realistic to take the credit 
quality of the financial institution for granted and to be highly superior to that of a general 
counterparty, no matter how prestigious or important the financial institution. 

Bilateral risk is also mentioned in the credit risk measurement space by the Basel II 
documentation. Annex IV, 2/A: "Unlike a firms exposure to credit risk through a loan, 
where the exposure to credit risk is unilateral and only the lending bank faces the risk of 
loss, the counterparty credit risk creates a bilateral risk of loss: the market value of the 
transaction can be positive or negative to either counterparty to the transaction." 

Basel II is more concerned with Risk Measurement than pricing. For an analysis of 
Counterparty risk in the risk measurement space we refer for example to De Frisco and Rosen 
(2005) [20j, who consider modeling of stochastic credit exposures for derivatives portfolios. 
However, also in the valuation space, bilateral features are quite relevant and often can be 
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responsible for seemingly paradoxical statementsjj For example, Citigroup in its press release 
on the first quarter revenues of 2009 reported a positive mark to market due to its worsened 
credit quality: "Revenues also included [...] a net 2.5USD billion positive CVA on derivative 
positions, excluding monolines, mainly due to the widening of Citis CDS spreads". In this 
paper we explain precisely how such a situation may origin. 

This paper first generalizes to the bilateral symmetric case the asymmetric unilateral 
setting proposed in [10], [8], [H] and [13], where only the counterparty is subject to default 
risk, while the investor is assumed to be default free. We provide a general formula that 
gives the bilateral risk credit valuation adjustment (BR-CVA) for portfolios exchanged be- 
tween a default risky investor and a default risky counterparty. Such formula shows that 
the adjustment to the investor is the difference between two discounted options terms, a 
discounted call option in scenarios of early default time of the counterparty minus a dis- 
counted put option in scenarios of early default times of the investor, both options being on 
the residual net present value of the portfolio at the relevant default times and having zero 
strike. The BR-CVA seen from the point of view of the counterparty is exactly the opposite. 
We allow for correlation between default of the investor, default of the counterparty and 
underlying portfolio risk factors, and for volatilities and dynamics in the credit spreads and 
in the underlying portfolio, all arbitrage free. 

Bilateral Counterparty adjustment applied to CDS. 

We then specialize our analysis to Credit Default Swaps (CDS) as underlying portfolio, 
generalizing the work of Brigo and Chourdakis [10] who deal with unilateral and asymmetric 
counterparty risk for these contracts. Featuring a CDS as underlying, a third default time 
enters the picture, namely the default time for the reference credit of the CDS. We therefore 
assume that all three entities are subject to default risk and that the default events of 
investor, counterparty and reference credit are correlated. We then propose a numerical 
methodology to evaluate the resulting BR-CVA formula. We investigate the impact of both 
credit spread volatility and default correlation on the credit valuation adjustment. Most of 
previous approaches on CDS counterparty risk only focus on unilateral counterparty risk, 
moreover ignoring the effect of volatility on the adjustment and mainly focusing on default 
correlation. 

The few earlier works on CDS with counterparty risk include Leung and Kwok (2005) [24j 
who, building on Collin-Dufresne et al. (2002) [16], model default intensities as deterministic 
constants with default indicators of other names as feeds. The exponential triggers of the 
default times are taken to be independent and default correlation results from the cross feeds, 
although again there is no explicit modeling of credit spread volatility. Hull and White [23] 
study the counterparty risk problem by resorting to barrier correlated models. Walker (2005) 
|30j models CDS counterparty risk using transitions rates as natural means to represent 
contagion, but again ignores credit spread volatility. Hille et al. (2006) [22] concentrate on 
credit risk measurement for CDS rather than precise valuation under counterparty risk. 

The need for explicitly modeling credit spread volatility is even more pronounced if the 
underlying reference contract is itself a CDS, as the credit valuation adjustment would involve 
CDS options and it is very undesirable to model options without volatihty in the underlying 



We are grateful to Dan Rosen for first signaling this issue to us during a conference in June 2009 
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asset. This has been clearly shown in [TO] for CDS, and in [T3] for interest rates payoffs, and 
again in [8] for commodities payoffs. Indeed, most credit models in the industry, especially 
when applied to Collateralized Debt Obligations or k-th to default baskets, model default 
correlation but ignore credit spread volatility. Credit spreads are typically assumed to be 
deterministic and a copula is postulated on the exponential triggers of the default times to 
model default correlation. This is the opposite of what used to happen with counterparty 
risk for interest rate underlyings, for example in Sorensen and Bollier (1994) [28] or Brigo 
and Masetti [TT] (in Pykhtin (2006) [27J), where correlation was ignored and volatility was 
modeled instead. Brigo and Chourdakis (2008) [TU] rectify this when addressing CDS's, but 
only deal with unilateral and asymmetric counterparty risk. 

Summarizing, we can put this paper's contribution into context with respect to its analo- 
gous earlier versions for unilateral counterparty risk or for other asset classes through Table [1} 

Other recent works include Crepey et al (2009) [19j who model wrong way risk for CDS 
with counterparty risk using a Markov chain copula model, and Blanchet-Scalliet and Patras 

(2008) [3], who resort to older Merton type models. A structural model with jumps has been 
introduced by Lipton and Sepp (2009) [25]. In these models assumptions on credit spread 
volatility are most time implicit. 

Here instead we introduce stochastic intensity models for all three names, including the 
investor. We do not introduce correlation between the three intensity processes, but model 
default correlation through a trivariate copula function on the exponential triggers of the 
default times. This is because typically spread correlation has a much lower impact on 
dependence of default times than default correlation. The correlation structure underlying 
the trivariate copula function may be estimated by the implied correlation in the quoted 
indices tranches markets such as i-Traxx and CDX or calculated from the asset correlation 
of the names. We leave this problem for further research. We present a preliminary numerical 
investigation to highlight the impact of dynamics parameters on bilateral CVA, followed by 
a case study of bilateral risk based on British Airways, Lehman Brothers and Royal Dutch 
Shell. We confirm findings in [TOj showing that both default correlation and credit spread 
volatilities have a relevant and structured impact on the adjustment. 

We specify that we do not consider specific collateral clauses or guarantees in the present 
work. We assume we are dealing with counterparty risk for an over the counter CDS trans- 
action where there is no periodic margining or collateral posting. The aim of this paper is 
analyzing the fine structure of counterparty risk adjustments with respect to finely tuned 
market dynamics, including the fundamental credit spread volatility ignored in other ap- 
proaches, and wrong way risk. We plan to address collateral provisions inclusion in future 
work. We already addressed netting in the interest rate context in and [TT]. The impact 
of credit triggers for the counterparty on CVA are analyzed in Yi (2009) [2^. Assefa et al 

(2009) [Ij analyze the modeling of coUateralization and margining in CVA calculations. 

Paper structure. 

The rest of the paper is organized as follows. Section |2] gives the general BR-CVA formula 
and refers to the appendix for the mathematical proof. Section |3] defines the underlying 
framework needed for applying the above methodology to credit default swaps, including 
the stochastic intensity model and the trivariate copula function for correlating defaults. 
Section |4] develops a numerical method to implement the formula in the case when the un- 
derlying contract is a CDS and gives the pseudo code of the algorithm used to run numerical 
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Table 1: Part of earlier analogous literature on CVA valuation with respect to inclusion 
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simulations. Section |5] present numerical experiments on different default correlation sce- 
narios and changes in credit spread volatility. Section |6] presents a specific application for 
computing the mark-to-market value of a CDS contract agreed between British Airways, 
Lehman Brothers and Royal Shell. Section [7| concludes the paper. 

2 Arbitrage- free valuation of bilateral counterparty risk 

We refer to the two names involved in the financial contract and subject to default risk as 

investor — name "0" 
counterparty name "2" 

In cases where the portfolio exchanged by the two parties is also a default sensitive 
instrument, we introduce a third name referring to the reference credit of that portfolio 

reference credit — > name "1" 

If the portfolio is not default sensitive then name "1" can be removed. 

We follow [6] and denote by tq, (ti) and T2 respectively the default times of the investor, 
(reference credit) and counterparty. We place ourselves in a probability space (fi, S, St,Q)- 
The filtration St models the flow of information of the whole market, including credit and 
Q is the risk neutral measure. This space is endowed also with a right-continuous and 
complete sub-filtration 3^t representing all the observable market quantities but the default 
event, thus 3't C Sf := 9^^ V JC*. Here, ^Kt = a({ro < u}{\/{ti < u}) V {t2 < u} : u < t) is the 
right-continuous filtration generated by the default events, either of the investor or of his 
counterparty (and of the reference credit if the underlying portfolio is credit sensitive). We 
also introduce the notion of stopped filtration. If r is an 3"^ stopping time, then the stopped 
filtration 5"^ is defined as 



S'r = a{Jt U {t < r},t > 0) 



(2.1) 
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If T is an St stopping time, then the stopped filtration Sr is defined as 

Sr = (T(StUO<T},t>0) (2.2) 

Let us call T the final maturity of the payoff which we need to evaluate and let us define 
the stopping time 

T = min{To,T2} (2.3) 

If T > T, there is neither default of the investor, nor of his counterparty during the hfe 

of the contract and they both can fulfill the agreements of the contract. On the contrary, if 
T < T then either the investor or his counterparty (or both) defaults. At r, the Net Present 
Value (NPV) of the residual payoff until maturity is computed. We then distinguish two 
cases: 

• T = T2. If the NPV is negative (respectively positive) for the investor (defaulted 
counterparty), it is completely paid (received) by the investor (defaulted counterparty) 
itself. If the NPV is positive (negative) for the investor (counterparty) , only a recovery 
fraction Rbc,2 of the NPV is exchanged. 

• T = Tq. If the NPV is positive (respectively negative) for the defaulted investor (coun- 
terparty), it is completely received (paid) by the defaulted investor (counterparty) 
itself. If the NPV is negative (positive) for the defaulted investor (counterparty), only 
a recovery fraction Rec.o of the NPV is exchanged. 

Let us define the following (mutually exclusive and exhaustive) events ordering the default 
times 



A = {to<T2<T} E^{T<To<T2} 
B = {To<r<T2} F = {r<T2<To} 

C = {T2<ro<r} 

D = {T2<r<To} (2.4) 

Let us call n^(t, T) the discounted payoff of a generic defaultable claim at t and Cashflows(m, s) 
the net cash fiows of the claim without default between time u and time s, discounted back 
at M, all payoffs seen from the point of view of the investor. Then, we have NPV(Ti) = 
E^.{CASHFLows(rj, T)}, i = 0,2. Let us denote by D{t,T) the price of a zero coupon bond 
with maturity T. We have 

Yi^it^T^ — l£;ijpCASHFLOWs(t, T") 

[CASHFLOWs(i,T2) +L>(i,T2) (Rbc.2(NPV(t2))+ - ( "NP V(t2) ) +) ] 

+1aub [CASHFi.ows(i, To) + D{t, To) ((NPV(to))+ - Rbco (-NPV(to))+)] 

(2.5) 

This last expression is the general payoff under bilateral counterparty default risk. Indeed, if 
there is no early default, this expression reduces to risk neutral valuation of the payoff (first 
term on the right hand side) . In case of early default of the counterparty, the payments due 
before default occurs are received (second term), and then if the residual net present value 



D. Brigo and A. Capponi, Arbitrage free bilateral counterparty risk valuation and application to CDS 7 



is positive only the recovery value of the counterparty Rec,2 is received (third term), whereas 
if it is negative it is paid in full by the investor (fourth term). In case of early default of 
the investor, the payments due before default occurs are received (fifth term), and then if 
the residual net present value is positive it is paid in full by the counterparty to the investor 
(sixth term), whereas if it is negative only the recovery value of the investor Rec,o is paid to 
the counterparty (seventh term). 

Let us denote by 11 (t, T) the discounted payoff for an equivalent claim with a default-free 
counterparty, i.e. n(t,T) = CASHFLows(t, T). We then have the following 

Proposition 2.1 ([6J). (General bilateral counterparty risk pricing formula). At 

valuation time t, and conditional on the event {t > t}, the price of the payoff under bilateral 
counterparty risk is 

Ei{n^(t,T)} = Et{U{t,T)} 

+Et {LGDo ■ Iaub ■ D{t,To) ■ [-NPV(ro)] + } 

-Et {LGD2 ■ IcuD ■ D{t,r2) ■ [NPV(r2)] + } (2.6) 

where LGD = 1 — Rec is the Loss Given Default and the recovery fraction Rec can he 
stochastic and possibly correlated with the default indicator process. It is clear that the value 
of a defaultable claim is the value of the corresponding default-free claim plus a long position 
in a put option (with zero strike) on the residual NPV giving nonzero contribution only in 
scenarios where the investor is the earliest to default (and does so before final maturity) 
plus a short position in a call option (with zero strike) on the residual NPV giving nonzero 
contribution in scenarios where the counterparty is the earliest to default ( and does so before 
final maturity), 



Proposition 2.1 is stated in ^ without a proof. Here, we provide a mathematical proof in 
Appendix|A| The adjustment is called bilateral counterparty risk credit valuation adjustment 
(BR-CVA) and it may be either positive or negative depending on whether the counterparty 
is more or less likely to default than the investor and on the volatilities and correlation. The 
mathematical expression is given by 

BR-CVA(t,T,LGDo,i,2) = {LGD2 ■ Icud ■ D{t,T2) ■ [NPV(r2)] + } 

-Et {LGDo ■ Iavjb ■ D{t, tq) ■ [-NPV(ro)]+} (2.7) 



where the right hand side in Eq. (2.7) depends on T through the events A,B,C,D and 
LGD012 = (LGDo, LGDi, LGD2) is shorthand notation to denote the dependence on the loss 
given defaults of the three names. 

Remark 2.2. (Symmetry vs Asymmetry). With respect to earlier results on counter- 



party risk valuation, Equation (2.7) has the great advantage of being symmetric. This is to 
say that if "2" were to compute counterparty risk of her position towards "1", she would 
find exactly — BR-CVA(t, T, LGDo,i,2)- However, if each party computed the adjustment by 
assuming itself to be default-free and considering only the default of the other party, then the 
adjustment calculated by "0" would be 

Et {LGD2 ■ lr,<T ■ D{t,T2) ■ [NPV(r2)] + } 

whereas the adjustment calculated by "2" would be 

Et {LGDo • lr,<T ■ D{t,To) ■ [-NPV(ro)] + } 
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and they would not be one the opposite of the other. This means that only in the first case 
the two parties agree on the value of the counterparty risk adjustment. 

Remark 2.3. (Change in sign). Earlier results on counterparty risk valuation, concerned 
with a default-free investor and asymmetric, would find an adjustment to be subtracted that 
is always positive. However, in our symmetric case even if the initial adjustment is positive 
due to 

Et {LGD2 • IcuD ■ D{t,r2) ■ [NPV(r2)] + } > {LGDo ■ Uub ■ D{t,To) ■ [-NPV(ro)] + } 

the situation may change in time, to the point that the two terms may cancel or that the 
adjustment may change sign as the credit quality of "0" deteriorates and that of "2" improves, 
so that the inequality changes direction. 

Remark 2.4. (Worsening of credit quality and positive mark to market). // the 



Investor marks to market her position at a later time using Formula (2.6), we can see that 
the term in Lgdq increases, ceteris paribus, if the credit quality of "0" worsens. Indeed, if we 
increase the credit spreads of the investor, now tq < T2 will happen more often, giving more 
weight to the term in Lgdq. This is at the basis of statements like the above one of Citigroup. 



3 Application to Credit Default Swaps 



In this section we use the formula developed in Section |2] to evaluate the BR-CVA in credit 



default swap contracts (CDS). Subsection 3.1 recalls the general formula for CDS evaluation. 



Subsection |3.2| introduces the copula models used to correlate the default events. Subsection 
3.3| recalls the CIR model used for the stochastic intensity of the three names. Subsection 
3.4| applies the general BR-CVA formula to calculate the adjustment for CDS contracts. 
We restrict our attention to CDS contract without an upfront trading. However, the pro- 
posed methodology is perfectly applicable to the case when the CDS trades with an upfront 
premium, which has become the case after the big bang protocol. 



3.1 CDS Payoff 

We assume deterministic interest rates, which leads to independence between ri and D{0, t), 
and deterministic recovery rates. Our results hold also true for the case of stochastic rates 
independent of default times. The receiver CDS valuation, for a CDS selling protection 
LGDi at time for default of the reference entity between times and Tf, in exchange of a 
periodic premium rate 5*1 is given by 



CDS,,fc(0,5i,LGDi) 



^1 



D{0,t){t-T^^t)-i 



+ J2 «i^(o,T, 

Tt 



>Ti) 



>t) 



i=a+l 

-t-LGDi 



D{0,t)dQ{Ti > t) 



where 7(t) is the first payment period Tj following time t. Let us denote by 

NPV(T,-,Tfe) := CDS„,6(T,-,^,LGDi) 



(3.1) 



(3.2) 
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the residual NPV of a receiver CDS between Ta and evaluated at time Tj, with Ta < Tj < 
Tfj. Eq. (3.2) can be written similarly to Eq. (3.1), except that now evaluation occurs at 
time Tj and has to be conditioned on the information set available to the market at Tj . This 
leads to 



CDS„,„(T,-,5i,LGDi) 



lri>r,CDSa,fe(Tj, Si,LGDi 



Lri>r, 



Si 



max{Ta,Tj} 



DiTj,t)it-T,^t)_, 



+ a,D{Tj,Ti)Q{Ti>T,\'5T, 

i=max{a 

+LGDi D{T,,t)dQ{ri>t\9T,] 



>t\S 



(3.3) 



For conversion of these running CDS into upfront ones, following the so called Big Bang 
protocol by ISDA, see for example Beumee, Brigo, Schiemert and Stoyle (2009) Our 
reasoning still applies with obvious modifications to the upfront CDS contract. 



3.2 Default Correlation 

We consider a reduced form model that is stochastic in the default intensity for the investor, 
counterparty and CDS reference credit. The default correlation between the three names is 
defined through a dependence structure on the exponential random variables characterizing 
the default times of the three names. Such dependence structure is modeled using a trivariate 
copula function. Let us denote by Xi{t) and Aj(t) = \i{s)ds respectively the default 
intensity and cumulated intensity of name i evaluated at time t. We recall that i = refers 
to the investor, i = 1 refers to the reference credit and i = 2 to the counterparty. We assume 
Aj to be independent of Xj for i j, and assume each of them to be strictly positive almost 
everywhere, thus implying that Aj is invertible. We stress the fact that independence of A's 
across names does not mean that the default event of one name does not change the default 
probability or intensity of other names, as discussed in [10], Section 4. We place ourselves 
in a Cox process setting, where 

r, = A-i(e.), z = 0,l,2 (3.4) 

with ^1 and .^2 being standard (unit-mean) exponential random variables whose associated 
uniforms 

Ui = l- exp{-^i} (3.5) 
are correlated through a Gaussian trivariate copula function 

Cr{uo, Ml, M2) = Q(f/o < Mo, f/l < Ml, f/2 < U2) (3.6) 

with R = [rjj]jj=o,i,2 being the correlation matrix parameterizing the Gaussian copula. No- 
tice that a trivariate Gaussian copula implies bivariate Gaussian marginal copulas. We show 
it for the case of the bivariate copula connecting the reference credit and the counterparty, 
but the same argument applies to the other bivariate copulas. Let {Xi, X2, X^) be a stan- 
dard Gaussian vector with correlation matrix R and let $jt be the distribution function of 
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a multivariate Gaussian random variable with correlation matrix R. For any pair of indices 
^ 7^ J; < j < 2, we denote by Rij the 2 ■ 2 submatrix formed by the intersection of row i 
and row j with column i and column j. We next state, without proof, an obvious result as 
a Lemma, which will be later used. 

Lemma 3.1. A trivariate Gaussian copula with correlation matrix R induces marginal bi- 
variate Gaussian copulas. 

We denote by Cij{ui,Uj) the bivariate copula associated to Rij. 
3.3 CIR stochastic intensity model 

We assume the following stochastic intensity model [7], [9] for the three names 

\j{t) = y,it)+^,{t;(3,), t>0,j = 0,1,2 (3.7) 

where ip is a deterministic function, depending on the parameter vector f3 (which includes 
Uo), that is integrable on closed intervals. We assume each yj to be a Cox IngersoU Ross 
(CIR) process [12] given by 

dyj{t) = K,{ii^ - yj{t))dt + uj^/^)dZ,{t) + JMujdMj{t),j = 0, 1, 2 (3.8) 

where J's are i.i.d. positive jump sizes that are exponentially distributed with mean Xj and 
Mj are Poisson processes with intensity rrij measuring the arrival of jumps in the intensity 
Xj. The parameter vectors are f3j = {nj, fij,h'j,yj{0),Xj,fnj) with each vector component 
being a positive deterministic constant. We relax the condition of inaccessibility of the origin 
2fvj/ij > Uj so that we do not limit the CDS implied volatility generated by the model. We 
assume the Zj^s to be standard Brownian motion processes under the risk neutral measure. 
We define the following integrated quantities which will be extensively used in the remainder 
of the paper 

Aj{t) = [ Xj{s)ds, Yj{t) = [ yj{s)ds, ^j{t;l3j)= [ ipjis; l3j)ds (3.9) 
Jo Jo Jo 

In this paper, we focus on intensities without jumps, i.e. mj = 0. Brigo and El-Bachir 

[9] consider in detail the tractable model with jumps, and this extension will be applied in 

future work. 



3.4 Bilateral risk credit valuation adjustment for receiver CDS 

We next proceed to the valuation of the BR-CVA adjustment for the case of CDS payoff 
given by Eq. (3.1). We state the result as a Proposition. 

Proposition 3.2. The BR-CVA at time t for a receiver GDS contract (protection seller) 
running from time Ta to time Ti, with premium S is given by 

BR-CVA-CDS,,b(t,5,LGDo,i,2) = LGD2 -Et {Icud ■ D{t,T2) ■ [lr,>r, CDSa,bir2, S, LGD,)]''^ 

-LGDo • Et {UuB ■ D{t,To) ■ [-lr,>roCDSa,b{ro,S,LGDi)]^} 

(3.10) 
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Proof. We have 

Et {icuD ■ D{t,T2) ■ [NPV(r2)] + } = Et {Icud ■ D{t,T2) ■ [CDS,,fc(r2, 5, LGDi)] + } 

IcuD ■ D{t, T2) ■ [lri>T2CDSa,6(r2, S, LGDi 



Et 



(3.11) 



where the first equahty in Eq. (3.11) follows by definition, while the last equality follows 
from Eq. (3.3). Similarly, we have 

{Iaub ■ D{t,T^) ■ [-NPV(ro)] + } = [Ia^b ■ D{t,ro) ■ [-CDS,,b(ro, 5, LGDi)] + } 

= Et{lAuB-D{t,ro) ■ [-1. 



'-Tl>TO 



CDSa,fe(ro,5,LGD] 



The proof follows using the expression of BR-CVA which is given by Eq. (2.7). 



(3.12) 

□ 



To summarize, in order to compute the counterparty risk adjustment, we determine the 
value of the CDS contract on the reference credit "1" at the point in time T2 at which the 
counterparty "2" defaults. The reference name "1" has survived this point and there is a 
bivariate copula Ci 2 which connects the default times of the reference credit and of the 
counterparty "2". Similarly, in order to compute the investor risk adjustment, we determine 
the value of the CDS contract on the reference credit "1" at the point in time tq at which 
the counterparty "0" defaults. The reference name "1" has survived this point and there is 
a bivariate copula Cq 1 which connects the default times of the reference credit and of the 



in order to compute (3.11) and (3.12) are 



investor "0". It is clear from Eq. (3.3) and Eq. (2.7) that the only terms we need to know 



lcuDlTi>T2Q(ri > t\Sr2; 



(3.13) 



and 



lAUBlri>To 



>t\9J (3.14) 
In the next section, we generalize the numerical method proposed in |10j to calculate 



quantities (3.13) and (3.14) for the bilateral counterparty case. 



4 Monte-Carlo Evaluation of the BR-CVA adjustment 

We propose a numerical method based on Monte-Carlo simulations to calculate the BR- 
CVA for the case of CDS contracts. Subsection |4.1| specifies the simulation method used to 



generate the sample paths of the CIR process. Subsection 4.2 gives a method to calculate 



CVA in (3.10) 



(3.13), while Subsection 4.3 gives the complete numerical algorithm for calculating the BR- 



4.1 Simulation of CIR process 

We use the well known fact [12] that the distribution of y(t) given y{u), for some m < t is, 
up to a scale factor, a noncentral chi-square distribution. More precisely, the transition law 
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of y{t) given y{u) can be expressed as 



yit) 



Ak 



where 



d 



(4.1) 



(4.2) 



and Xui"^) denotes a non-central chi-square random variable with u degrees of freedom and 
non centrality parameter v. In this way, if we know 1/(0), we can simulate the process y{t) 
exactly on a discrete time grid by sampling from the non-central chi-square distribution. 

4.2 Calculation of Survival Probability 

We state the result in the form of a proposition. Such result will be used in Subsection 4^ 
to develop a numerical algorithm for computing the BR-CVA. Let us define 



exp(-Ai(rj)) 



(4.3) 



and denote by Fx-(^t) the cumulative distribution function of the cumulative (shifted) intensity 
of the CIR process associated to name i, which can be retrieved inverting the characteristic 
function of the integrated CIR process. We have 

Proposition 4.1. 

lcuDlri>r2Q(n > t\Sr2) = 



lr2<Tlr2<ro + '^T2<t^Tx>T2 



^Ai(t)-Ai(T2)(- log(l - Ui) - Al(r2))c?Ci|o,2(tii; U2) 

(44) 



where 



A 



{t<T2< Ti} 

9Ci,2(«l,M2) 



du2 



Ci|o,2(mi; U2) 



dC(Uo^2,ui,U2) 



U2 = U2 



du2 



8C1, 2(^1, 2,^2) 



U2 = U2 



du2 



+ 



dC{Uo,2,Ux,2,U2) 



U2 = U2 



du2 



U2=U2 











dCo,2{Ul\2,U2) 

du2 


9Cl,2(f^l,2,M2) 

du2 


1 dC{Uo,2,Ui,2,U2) 

du2 






U2=U2 


U2=U2 


U2 = U2 



(4.5) 



Similarly, 
Proposition 4.2. 

lAUBlri>ToQ(Tl > ^ISro) = 
lro<Tlro<T2 1_B + lro<tlri>ro 



"Ai(i)-Ai(To)l 



■log(l - Ml) - Ai(ro))rfCi|2,o(Mi; Uo) 



(4.6) 
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where 



B = {t<ro< Ti} 



Ci|2,o(mi; Uq) 



dCo^i{uo,ui) 
duo 



dC{uo,ui,U2,o) 



uo=Uo 



duo 



uo=Uo 



dCo^i{uo,Uifi) 
duo 



uo=Uo 



dC{uo,Ul,o,U2.o) 

duo 



uo=Uo 











dCo,2{uo,U2.o) 

duo 


dCo,i(uo,Uifi) 
duo 


dCiuo,Ui,o,U2,o) 

duo 






uo=Uo 


uo=Uo 


uo=Uo 



(4.7) 



The proofs of the propositions are reported in Appendix [B] 

4.3 The numerical BR-CVA adjustment algorithm 

We give the pseudo code of the numerical algorithm used to calculate the BR-CVA for the 
CDS payer and receiver. In the pseudo-code below, the variable represents the time 
elapsing between payment period and measured in years, the variable A represents 



the fineness of the grid used to evaluate the integral of the survival probability in Eq. (4.4) 
and the variable 6 represents the fineness of the time grid used to evaluate the integral in 
Eq. (3.3). The variable Xmax represents the maximum x value for which the cumulative 
distribution function is implied from the characteristic function. The expression CDF(xfc) 
corresponds to a subroutine call which calculates the cumulative distribution function of the 
integrated CIR process at Xk- This can be done by inversion of the characteristic function 
of the integrated CIR process using Fourier transform methods such as in [18]. The inputs 
to the main procedure Calculate Adjustment are the number of Monte-Carlo runs 
and the market quote 5*1 of the 5 year CDS spread of the reference entity. 



Algorithm 1 [BR-CVA_R, BR-CVA_P] 
for z = 1 : A^ do 



Calculate Adjustment (N, S) 



Generate tq, ti, and T2 using Eq. (3.5) and Eq. (3.6). 
if T2 < Tq and T2 < Tf, then 
if Ti > T2 then 

[BR-CVA_R_2, BR-CVA_P_2] = CDSAdjust(T^(,2), S, LGDi, 2) 
CUM_BR-CVA_R = CUM_BR-CVA_R + LGDa- BR-CVA_R_2 
CUM_BR-CVA_P = CUM:BR-CVA_P + LGDa- BR-CVA_P_2 
end if 
end if 

if To < T2 and tq < then 
if Ti > To then 

[BR-CVA_R_0, BR-CVA_P_0] = CDSAdjust(T^(,„), S, LGDi, 0) 



CUM_BR-CVA_R 
CUM_BR-CVA_P ■ 



CUM_BR-CVA_R 
CUM_BR-CVA_P 
end if 
end if 
end for 

BR-CVA_R = CUM_BR-CVA_R / A^ 
BR-CVA_P = CUM_BR-CVA_P / A^ 



7{ to); 

LGDo- BR-CVA_R_0 
LGDo- BR-CVA_P_0 
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Algorithm 2 [CDSr,CDSp] = CDSAdjust(T,, 5, LGDi, index) 



Terrrii = Term2 = Term^ = 

t start = Kiax(Ta,Tj) 

Qprev = ComputeProb (tstart, Tj, index) 
for t = t start + 5 : 5 -.Tbdo 

Qcurr = ComputeP Toh(t,Tj, index) 

Termi = Ternii + D{Tj,t - 6)(t - 6 - T^{t~S)~i) {Qcurr - Qprev) 

Terms = Term^ + D{Tj,t- 6) {Qcurr - Qprev) 
Qprev Qcurr 

end for 

for ti = tstart + ai'. ai'.Tb do 

Qcurr = ComputeP Toh{ti,Tj, index) 

Term2 = Term2 + (Xi ■ D{Tj, ti) ■ Qcurr 
end for 

CDSvai = S ■ {Term2 — Ternii) + LGDi ■ Terms 
if index == 2 then 

CDSr = D{t, Tj) ■ max{CDS^ah 0) 

CDSp = D{t, Tj) ■ max{-CDSyah 0) 
end if 

if index == then 

CDSr = D{t, Tj) ■ max{-CDS,ai, 0) 

CDSp = D{t,Tj) ■ max{CDS^ai,0) 
end if 



Algorithm 3 Qi = CoMPUTEPROB(t, T,-, mdea;) 



Uindex 1 6Xp{ Y^index{Tj) ^ index{Tj j /3index)} 

U, = l- exp{-Y,{T,) - 
for Xfe = : A : x^ax do 
Pk = CDF(xfc) 

Uk = l- exp{-XA: - ^l(t)} 

if index == 2 then 

Compute fk = Ci|o,2(tifc; U2) using Eq. 1^^^ 
else 

Compute fk = Ci\2,o{uk; Uq) using Eq. ( |47| ) 
end if 
end for 

Qi = 12{uk,pkJk)--nk>UiPk{fk+l " fk) 



5 Numerical Results 

We consider an investor (name "0") trading a five- years CDS contract on a reference name 
(name "1") with a counterparty (name "2"). Both the investor and the counterparty are 
subject to default risk. We experiment on different levels of credit risk and credit risk 
volatility of the three names, which are specified by the parameters of the CIR processes in 
Table [2j We recall that the survival probabilities associated with a CIR intensity process 
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Credit Risk Levels 


1/(0) 






Credit Risk volatilities 


V 


low 


0.00001 


0.9 


0.0001 


low 


0.01 


middle 


0.01 


0.80 


0.02 


middle 


0.2 


high 


0.03 


0.50 


0.05 


high 


0.5 



Table 2: The credit risk levels and credit risk volatilities parameterizing the CIR processes 



Maturity 


Low Risk 


Middle Risk 


High risk 


ly 





92 


234 


2y 





104 


244 


3y 





112 


248 


4y 


1 


117 


250 


5y 


1 


120 


251 


6y 


1 


122 


252 


7y 


1 


124 


253 


8y 


1 


125 


253 


9y 


1 


126 


254 


lOy 


1 


127 


254 



Table 3: Break-even spreads in basis points generated using the parameters of the CIR 
processes in Table [2j The first column is generated using low credit risk and credit risk 
volatility. The second column is generated using middle credit risk and credit risk volatility. 
The third column is generated using high credit risk and credit risk volatility. 



are given by 



[r,>t) := E[e-^^W] 

= P^^^(0,t,A) (5.1) 



where P^-^^(0, t, A) is the price at time of a zero coupon bond maturing at time t under a 
stochastic interest rate dynamics given by the CIR process ^2], with A = {yi{0), Ki, fii,h'i) 
being the vector of CIR parameters, i = 0,1, 2. 



We report in Table 3 the break-even spreads zeroing (3.1) in S*, with survival probabilities 



given by Eq. (5.1) and CIR parameters l3iow, (^middle, (^high obtained from Table |2| 

The evaluation time t and the starting time of the CDS contract are both set to 
zero. The end time of the contract is set to five years. It is assumed that payments 
are exchanged every three months. The loss given defaults of the low, middle and high 
risk entity are respectively set to LGD^o^ = 0.6, hGDmiddie = 0.65, LGDhigh = 0.7. We 
assume that the spreads in Table |3] are the spreads quoted in the markets for the three 
names under consideration. We recover the integrated shift \E'(t; f3) which makes the model 
survival probabilities consistent with the market survival probabilities coming from Table 
|3] whenever we change the CIR parameters. In mathematical terms, for any i = 0, 1,2, we 
impose that 

Q{Ti>t)^odei ■■= E[e^^'W] 

= Qi^i > t)market (5-2) 
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The market survival probability for name "z" is bootstrapped from the market CDS quotes 
reported in Table [3] Such bootstrap procedure is performed assuming a piecewise linear 
hazard rate function. From the definition of the integrated process '^j{t](3j) given in Eq. 



(3.9), we can restate Eq. (5.2) as 

^,(t;A) = "^ogi—^ - 



market 



,0, , , ,5,, 



[n > t) 



market 



where the last equality in Eq. (5.3) follows from Eq. (5.1). We first study a case where 
the quoted market spreads of the investor has a low credit risk profile (Piow), the reference 
entity has high credit risk profile [Phigh), while the counterparty has middle credit risk 
iPmiddie)- We Vary the correlation between reference credit and counterparty as well as the 
credit risk volatility I'l of the reference credit. Since the focus is mostly on credit spreads 
volatility, we calculate the implied CDS volatility produced by the choice of parameters /3i = 
(l/o(l), /ii, Ki, i^i) for a hypothetical CDS option, maturing in one year and in case of exercise 
entering into a four year CDS contract on the underlying reference credit "1" . The objective 
of the experiments is to measure the impact of correlation and credit spreads volatility 
on the BR-CVA. The triple {x, y, z) represents the correlation of the trivariate Gaussian 
copula, with x denoting the correlation between the investor and reference credit, y denoting 
the correlation between the investor and the counterparty and z denoting the correlation 
between the reference credit and the counterparty. The values BR-CVA_P and BR-CVA_R 
are respectively the Monte-Carlo estimates of the CDS payer and receiver counterparty risk 
adjustments. The theoretical formula for payer risk adjustment is given by 

LGD2 ■ Et {icuD ■ D(t, T2) ■ [-NPV(r2)] + } 

-LGDo ■ Et {Iaub ■ D(t, tq) • [NPV(ro)] + } (5.4) 



which follows from adapting the formula given in Eq. (3.10), given for the receiver counter- 
party adjustment, to the case of the payer counterparty adjustment. Tables |4] and [5] report 
the results obtained. 

The scenarios considered in Table |4] and [5] assume an investor with an extremely low credit 
risk profile and thus are similar to those considered in [10], where the investor is assumed 
to be default-free. Therefore, it is not surprising that we find similar results. Similarly 
to them, we can see that the BR-CVA for the payer investor tends to be monotonically 
increasing with the correlation ri2, whereas the BR-CVA for the receiver investor appears to 
be monotonically decreasing. This is expected since for low and negative correlation values, 
defaults of the counterparty come with reductions in default risk of the reference entity, thus 
the receiver investor holds an option which is in the money, but at the counterparty default 
he only gets a fraction of it proportional to the recovery value of the counterparty. On the 
contrary, for high and positive correlation values, defaults of the counterparty come with 
increases in default risk of the reference entity, thus the payer investor holds an option which 
is in the money, but at the counterparty default he only gets a fraction of it proportional 
to the recovery value of the counterparty. While these patterns are intuitive, there is a 
particular pattern in correlation that is seemingly counterintuitive. 

Remark 5.1. (Decreasing wrong way risk with low credit spread volatility and 
high correlation). In general, we expect the BR-CVA for the payer investor to increase with 
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(roi,ro2,ri2) 


Vol. parameter vi 
CDS Impled vol 


0.01 
1.5% 


0.10 
15% 


0.20 
28% 


0.30 
37% 


0.40 
42% 


0.50 
42% 


(0, 0, -0.99) 


BR-CVA_P 
BR-CVA_R 


0.0(0.0) 
29.3(1.5) 


0.0(0.0) 
29.7(1.5) 


0.0(0.0) 
29.8(1.5) 


0.0(0.0) 
29.8(1.5) 


-0.0(0.0) 
29.8(1.5) 


-0.0(0.0) 
29.6(1.5) 


(0, 0, -0.90) 


BR-CVAJ* 
BR-CVAJl 


0.0(0.0) 
29.6(1.5) 


0.0(0.0) 
29.2(1.5) 


0.0(0.0) 
29.1(1.5) 


0.0(0.0) 
29.3(1.5) 


-0.0(0.0) 
29.6(1.5) 


-0.0(0.0) 
29.0(1.5) 


(0, 0, -0.60) 


BR-CVAJ* 
BR-CVAJl 


0.0(0.0) 
27.0(1.4) 


0.0(0.0) 
27.2(1.4) 


0.0(0.0) 
26.4(1.4) 


0.1(0.0) 
26.6(1.4) 


0.4(0.2) 
25.4(1.3) 


0.6(0.2) 
25.2(1.3) 


(0, 0, -0.20) 


BR-CVA_P 
BR-CVAJl 


0.0(0.0) 
10.3(0.7) 


0.4(0.1) 
10.4(0.7) 


1.3(0.2) 
10.9(0.7) 


2.0(0.3) 
12.1(0.7) 


2.6(0.4) 
12.8(0.8) 


3.2(0.5) 
13.0(0.8) 


(0, 0, 0) 


BR-CVAJ 

BR-CVA_R 


4.8(0.3) 

0.0(0.0) 


5.1(0.4) 

0.5(0.1) 


5.4(0.5) 

2.1(0.2) 


7.4(0.8) 
3.9(0.3) 


5.9(0.6) 
5.0(0.3) 


5.6(0.7) 

6.3(0.4) 


(0, 0, 0.20) 


BR-CVA_P 
BR-CVAJl 


25.9(1.5) 
0.0(0.0) 


25.1(1.5) 
0.0(0.0) 


23.2(1.5) 
0.2(0.0) 


21.1(1.6) 
0.5(0.1) 


15.9(1.3) 
0.7(0.1) 


12.8(1.2) 
1.1(0.1) 


(0, 0, 0.60) 


BR-CVAJ 

BR-CVAJl 


72.0(4.9) 

0.0(0.0) 


72.0(4.8) 

0.0(0.0) 


65.2(4.4) 

0.0(0.0) 


57.1(4.0) 

0.0(0.0) 


50.8(3.6) 

0.0(0.0) 


40.4(3.0) 

0.1(0.0) 


(0, 0, 0.90) 


BR-CVA_P 
BR-CVAJl 


68.4(6.1) 
0.0(0.0) 


73.8(6.3) 
0.0(0.0) 


69.2(5.8) 
0.0(0.0) 


65.2(5.4) 
0.0(0.0) 


61.5(5.0) 
0.0(0.0) 


62.2(4.9) 
0.0(0.0) 


(0, 0, 0.99) 


BR-CVAJ' 
BR-CVAJl 


13.2(2.8) 
0.0(0.0) 


28.4(4.0) 
0.0(0.0) 


39.4(4.4) 
0.1(0.0) 


51.9(4.9) 
0.1(0.0) 


54.8(5.0) 
0.1(0.0) 


67.6(5.4) 
0.4(0.3) 



Table 4: BR-CVA in basis points for the case when yo{0) = 0.0001, kq = 0.01, /Uq = 
0.001,1/0 = 0.01, yo{2) = 0.01, /t2 = 0.8, /ia = 0.02, i/g = 0.2. The name "1" has CIR 
parameters ?/o(l) = 0.03, /ti = 0.5, /xi = 0.05, while the credit spreads volatihty vi is varied 
across scenarios. The numbers within round brackets represent the Monte-Carlo standard 
error. The CDS contract on the reference credit has a five-years maturity. 
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(roi,ro2,ri2) 


Vol. parameter vi 
CDS Impled vol 


0.01 
1.5% 


0.10 
15% 


0.20 
28% 


0.30 
37% 


0.40 
42% 


0.50 
42% 


(0, 0, -0.99) 


BR-CVA_P 
BR-CVA_R 


0.0(0.0) 
28.8(1.4) 


0.0(0.0) 
29.2(1.4) 


0.0(0.0) 
28.3(1.4) 


0.0(0.0) 
28.0(1.4) 


0.0(0.0) 
29.3(1.4) 


0.0(0.0) 
28.9(1.4) 


(0, 0, -0.90) 


BR-CVAJ* 
BR-CVAJl 


0.0(0.0) 
28.9(1.4) 


0.0(0.0) 
29.0(1.4) 


0.0(0.0) 
28.2(1.4) 


0.0(0.0) 
28.7(1.4) 


0.1(0.0) 
28.8(1.4) 


0.1(0.0) 
28.9(1.4) 


(0, 0, -0.60) 


BR-CVAJ* 
BR-CVAJl 


0.0(0.0) 
26.8(1.3) 


0.0(0.0) 
26.5(1.3) 


0.1(0.1) 
25.7(1.3) 


0.1(0.0) 
25.1(1.3) 


0.8(0.2) 
24.9(1.3) 


0.2(0.1) 
25.0(1.3) 


(0, 0, -0.20) 


BR-CVA_P 
BR-CVAJl 


0.0(0.0) 
9.7(0.6) 


0.4(0.1) 
9.7(0.6) 


1.3(0.2) 
10.4(0.7) 


2.4(0.4) 
11.7(0.7) 


2.8(0.5) 
12.7(0.7) 


2.0(0.4) 
13.0(0.7) 


(0, 0, 0) 


BR-CVAJ 

BR-CVA_R 


4.8(0.2) 
0.0(0.0) 


5.3(0.4) 
0.5(0.1) 


6.0(0.5) 

2.1(0.2) 


6.6(0.7) 

3.8(0.3) 


5.2(0.7) 
5.0(0.3) 


5.1(0.7) 
6.0(0.3) 


(0, 0, 0.20) 


BR-CVA_P 
BR-CVAJl 


26.6(1.5) 
0.0(0.0) 


26.1(1.5) 
0.0(0.0) 


23.6(1.5) 
0.1(0.0) 


19.7(1.4) 
0.4(0.0) 


16.7(1.4) 
0.7(0.1) 


13.8(1.4) 
1.0(0.1) 


(0, 0, 0.60) 


BR-CVAJ 

BR-CVAJl 


76.4(4.8) 

0.0(0.0) 


74.1(4.6) 

0.0(0.0) 


68.1(4.4) 

0.0(0.0) 


60.8(4.0) 

0.0(0.0) 


52.2(3.7) 

0.1(0.0) 


42.3(3.1) 

0.1(0.0) 


(0, 0, 0.90) 


BR-CVA_P 
BR-CVAJl 


75.3(6.1) 
0.0(0.0) 


76.1(6.1) 
0.0(0.0) 


74.4(5.9) 
0.0(0.0) 


68.8(5.5) 
0.0(0.0) 


62.2(5.2) 
0.0(0.0) 


64.4(4.9) 
0.0(0.0) 


(0, 0, 0.99) 


BR-CVAJ' 
BR-CVAJl 


12.6(2.6) 
0.0(0.0) 


24.9(3.6) 
0.0(0.0) 


41.3(4.5) 
0.1(0.0) 


51.8(4.9) 
0.1(0.0) 


55.3(4.9) 
0.1(0.0) 


65.9(5.3) 
0.1(0.0) 



Table 5: BR-CVA in basis points for the case when |/o(0) = 0.0001, kq = 0.01, ^uq = 
0.001,1/0 = 0.01, yo{2) = 0.01,^2 = 0.8, /i2 = 0.02, 1/2 = 0.01. The name "1" has CIR 
parameters ?/o(l) = 0.03, ki = 0.5, /xi = 0.05, while the credit spreads volatility vi is varied 
across scenarios. The numbers within round brackets represent the Monte-Carlo standard 
error. The CDS contract on the reference credit has a five-years maturity. 
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Bilataral Risk Counterparty Value Adjustment-Payer Investor 

80 I < < < < < 1 1 r— 




Figure 1: Patterns of the BR-CVA of the CDS contract for payer and receiver investor and 
for low (0.1) and high (0.3) reference credit spreads volatihty Ui. We have that 1^2 = 0.01 
and z/q = 0.01 

correlation between the underlying reference credit and the counterparty default. However, 
for correlation increasing beyond relatively large values, the BR-CVA for the payer investor 
goes down instead, if the credit spreads volatility of the reference entity is small. Consider, 
for example, the case when the reference entity and the counterparty are 99% correlated, 
so that the exponential triggers and ^2 ore almost identical. If we consider the scenario 
where ui = 1/2 = 0.01, then the intensity of both processes are almost deterministic, with 
Ai > A2 having the name "1" high credit risk and name "2" middle credit risk. Thus, to 
a first approximation Ti = ^ is higher than T2 = and consequently the reference credit 
always defaults before the counterparty, resulting in no adjustment. Table\^has qualitatively 
a similar behavior. This is visible also in the right side of the graphs in the upper Figure [I| 
where the pattern goes down at the very end for correlations beyond 0. 9. 

This seemingly counterintuitive feature is improved if credit spread volatility becomes large 
(and in line with realistic CDS implied volatilities, see Brigo (2005)1^). This is evident from 
the last column of Table\^ and is due to the increased randomness in the default times, that 
now can cross each other in more scenarios. 

To further illustrate the transition from large (small) to small (large) adjustments for the 
receiver (payer) investor and illustrate the effect of correlation and credit spreads volatility, 
we display in Figure [T] the behavior of the bilateral risk adjustment for two different values 
of credit spreads volatility. 

Tables |6] and [7] report the BR-CVA adjustments for the Payer CDS under a set of five 
different riskiness scenarios. The assignments of credit risks in Table [2] to the three names 
determine the scenario in place. We have: 

• Scenario 1 {Base Scenario). The investor has low credit risk, the reference entity has 
high credit risk, and the counterparty has middle credit risk 

• Scenario 2 [Risky counterparty) . The investor has low credit risk, the reference entity 
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('^01,'^02,n2) 


Base Scenario 


Risky Counterparty 


Risky Investor 


Risky Ref 


Safe Ref 


(0, 0, 0) 


6.0(0.4) 


3.6(0.2) 


-0.8(0.0) 


-0.0(0.0) 


5.5(0.4) 


(0, 0, 0.1) 


15.1(0.8) 


12.6(0.5) 


-0.8(0.0) 


-0.0(0.0) 


13.2(1.0) 


(0, 0, 0.3) 


37.0(2.0) 


37.4(1.5) 


-0.8(0.1) 


0.1(0.0) 


34.6(2.1) 


(0, 0, 0.6) 


73.7(4.4) 


92.5(3.8) 


-0.3(0.5) 


0.7(0.1) 


67.3(4.3) 


(0, 0, 0.9) 


83.4(6.0) 


207.8(8.2) 


-0.8(0.0) 


1.5(0.3) 


75.7(5.8) 


(0, 0, 0.99) 


26.0(3.5) 


316.9(12.5) 


-0.8(0.0) 


1.8(0.5) 


22.6(3.4) 


(0, 0.1, 0) 


6.0(0.4) 


3.6(0.2) 


-0.8(0.0) 


-0.0(0.0) 


4.2(0.3) 


(0, 0.3, 0) 


5.8(0.4) 


3.6(0.2) 


-0.8(0.0) 


-0.0(0.0) 


3.7(0.3) 


(0, 0.6, 0) 


6.0(0.4) 


3.6(0.2) 


-0.8(0.1) 


-0.0(0.0) 


5.9(0.4) 


(0, 0.9, 0) 


0.6(4.5) 


3.5(0.2) 


-0.8(0.0) 


0.0(0.0) 


15.0(1.5) 


(0.1, 0, 0) 


6.0(0.4) 


3.6(0.2) 


-0.1(0.0) 


-0.0(0.0) 


6.5(0.4) 


(0.3, 0, 0) 


6.0(0.4) 


3.6(0.2) 


0.0(0.0) 


0.0(0.0) 


6.6(0.4) 


(0.6, 0, 0) 


6.0(0.4) 


3.6(0.2) 


0.0(0.0) 


0.0(0.0) 


6.5(0.4) 


(0.9, 0, 0) 


5.9(0.4) 


3.6(0.2) 


0.0(0.0) 


-0.1(0.0) 


6.6(0.4) 


(0.99, 0, 0) 


5.9(0.4) 


3.6(0.2) 


-1.8(0.1) 


-0.1(0.0) 


6.7(0.4) 



Table 6: BR-CVA under five different riskiness scenarios. The CIR volatilities are set to 
Uq = Ui = 1/2 = 0.1. The correlation triple has two non-zero entries 



('^01,'^02,?^20) 


Base Scenario 


Risky Counterparty 


Risky Investor 


Risky Ref 


Safe Ref 


(0.7, 0.4, 0.3) 


37.0(2.1) 


36.4(1.4) 


-0.0(0.0) 


-0.0(0.0) 


8.4(0.7) 


(0.7, 0.3, 0.4) 


49.0(2.8) 


51.4(2.1) 


0.1(0.1) 


0.2(0.0) 


25.8(1.7) 


(0.4, 0.7, 0.3) 


35.3(2.0) 


35.5(1.4) 


0.1(0.1) 


0.1(0.0) 


27.0(1.7) 


(0.4, 0.3, 0.7) 


79.8(5.0) 


117.4(4.9) 


0.6(0.5) 


0.6(0.1) 


74.0(4.9) 


(0.3, 0.4, 0.7) 


80.8(5.1) 


117.8(4.9) 


0.0(0.0) 


0.6(0.1) 


73.6(5.5) 


(0.3, 0.7, 0.4) 


45.8(2.7) 


50.0(2.0) 


0.0(0.0) 


0.4(0.0) 


50.9(3.3) 


(0.3, 0.3, 0.3) 


36.4(2.0) 


36.7(1.4) 


0.0(0.0) 


0.1(0.0) 


28.2(1.7) 


(0.4, 0.4, 0.4) 


47.9(2.8) 


50.8(2.1) 


0.0(0.0) 


0.2(0.0) 


39.0(2.4) 


(0.7, 0.7, 0.7) 


77.4(5.0) 


113.6(4.8) 


0.0(0.0) 


0.6(0.1) 


82.5(6.5) 



Table 7: BR-CVA under five different riskiness scenarios. The CIR volatilities are set to 
i>o = ui = 1/2 = 0.1. The correlation triple has all non-zero entries 



has middle credit risk, and the counterparty has high credit risk. 

• Scenario 3 {Risky investor). The counterparty has low credit risk, the reference entity 
has middle credit risk, and the investor has high credit risk. 

• Scenario 4 [Risky Ref). Both investor and counterparty have middle credit risk, while 
the reference entity has high credit risk. 

• Scenario 5 [Safe Ref). Both investor and counterparty have high credit risk, while the 
reference entity has low credit risk. 

Tables |6] and [7] clearly show the effect of the wrong way risk. For example, looking at the 
second column, we can see that as the correlation between counterparty and reference entity 
gets larger, the BR-CVA increases significantly. This is because (1) the counterparty is the 
riskiest name and (2) the high positive correlation makes the spread of the reference entity 
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larger at the counterparty default, thus the option on the residual NPV for the payer investor 
will be deep into the money and worth more. If the only correlation is between reference 
credit and counterparty, then the adjustment is the largest, while when the counterparty is 
less correlated with the reference entity, but also positively correlated with the investor (see 
second column of Table [?]), then the adjustment tends to decrease as the number of scenarios 
where the counterparty is the earliest to default is smaller. 

The first six rows of Table [6] show how much the BR-CVA adjustment is driven by corre- 
lation and credit spreads volatility. When the counterparty and reference entity are loosely 
correlated, then the adjustments in Scenario 1 and Scenario 2 are very similar, although the 
counterparty is riskier in Scenario 2. However, as the correlation increases, the BR-CVA 
adjustment in Scenario 2 becomes significantly larger than the corresponding adjustment in 
Scenario 1. Differently from Tables |4] and |5| however, we have that low credit risk volatility 
values has the opposite effect and amplify the adjustment; this is because for 99% correla- 
tion between reference credit and counterparty, we now have that the counterparty always 
defaults earlier being riskier. 

Consistently with the results in Table |4] and |5| no BR-CVA takes place in scenario 4. 
This is because the reference entity has the highest credit risk profile and since the credit 
risk volatilities of the three names are relatively low, the reference entity always defaults 
first, thus resulting in no adjustment taking place. 

6 Application to a market scenario 

We apply the methodology to calculate the mark-to-market price of a five-year CDS contract 
between British Airways (counterparty) and Lehman Brothers (investor) on the default of 
Royal Dutch Shell (reference credit). We consider two CDS contracts. In the first contract 
Lehman Brothers buys 5-year protection on Shell from British Airways on January 5, 2006. 
In the second contract, Lehman Brothers sells 5-year protection on Shell to British Airways 
on January 5, 2006. In both contracts, British Airways computes the mark-to-market value 
of the contract on May 1, 2008. We consider different correlation scenarios among the three 
names. The CDS quotes of the three names on those dates are reported in Tables [8] and |9} 



Maturity 


Royal Dutch Shell 


Lehman Brothers 


British Airways 


ly 


4 


6.8 


10 


2y 


5.8 


10.2 


23.2 


3y 


7.8 


14.4 


50.6 


4y 


10.1 


18.7 


80.2 


5y 


11.7 


23.2 


110 


6y 


15.8 


27.3.3 


129.5 


7y 


19.4 


30.5 


142.8 


8y 


20.5 


33.7 


153.6 


9y 


21 


36.5 


162.1 


lOy 


21.4 


38.6 


168.8 



Table 8: Market spread quotes in basis points for Royal Dutch Shell, Lehman Brothers and 
British Airways on January 5, 2006. 
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Maturity 


Royal Dutch Shell 


Lehman Brothers 


British Airways 


1 V 


24 


203 


1 51 




24 fi 




230 




2fi 4 


1 66 75 


275 


4v 


28 5 


152.25 


305 




tJ\J 


145 


335 


6v 


32.1 


136.3 


342 


7y 


33.6 


130 


347 


8y 


35.1 


125.8 


350.6 


9y 


36.3 


122.6 


353.3 


lOy 


37.2 


120 


355.5 



Table 9: Market spread quotes in basis points for Royal Dutch Shell, Lehman Brothers and 
British Airways on May 1, 2008. 



Credit Risk Levels (2006) 


1/(0) 




fi 


V 


Lehman Brothers (name "0") 


0.0001 


0.036 


0.0432 


0.0553 


Royal Dutch Shell (name "1") 


0.0001 


0.0394 


0.0219 


0.0192 


British Airways (name "2") 


0.00002 


0.0266 


0.2582 


0.0003 



Table 10: The CIR parameters of Lehman Brothers, Royal Dutch Shell and British Airways 
calibrated to the market quotes of CDS on January 5, 2006. 



The calibrated parameters for the CIR process dynamics on those dates are reported 



in Table 10 and Table 11 respectively. The calibration is done assuming zero shift and 
inaccessibility of the origin. To have a perfect calibration we should add a shift, but in this 
case the quality of fit is relatively good also without shift and this keeps the model simple. 
The absolute calibration errors 1 5'^^'! — S^^s "^"^^ | for the different names at the different dates 
range from less than one basis point to a maximum of 20 basis points for BA quotes about 
355 basis points and of 23 for Lehman quotes about 188 basis points. These errors can be 
zeroed by introducing a shift, but for illustrating the CVA features the lack of shift does not 
compromise the analysis. 

The procedure used for mark-to- market value valuation is detailed next: 

(a) We take the CDS quotes of British Airways, Lehman Brothers and Royal Dutch Shell 
on January 5, 2006 and calibrate the parameters of the CIR processes associated to the 
three names assuming zero shift and inaccessibility of the origin. The results obtained 



Credit Risk Levels (2008) 


y(o) 


K 




V 


Lehman Brothers (name "0") 


0.6611 


7.8788 


0.0208 


0.5722 


Royal Dutch Shell (name "1") 


0.003 


0.1835 


0.0089 


0.0057 


British Airways (name "2") 


0.00001 


0.6773 


0.0782 


0.2242 



Table 11: The CIR parameters of Lehman Brothers, Royal Dutch Shell and British Airways 
calibrated to the market quotes of CDS on May 1, 2008. 
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from the calibration are the ones reported in Table 10 



(b) We calculate the value of the five year risk-adjusted CDS contract starting at Ta = 
January 5, 2006 and ending five years later at Tb = January 5, 2011 as 

C7D5f,(T,, Su LGDo,i,2) = CDSa,b{Ta, LGDi) - BR-CVA-CDS,,b(T,, 5i, LGDo,i,2) 

(6.1) 

where 5*1 = 120 bps is the five-year spread quote of Royal Dutch Shell at time Ta, 
CDSa^fc(Ta, 5*1, LGDi) is the value of the equivalent CDS contract which does not ac- 



count for counterparty risk given by Eq. (3.1) and the loss given default of the three 



names are taken from a market provider and equal to 0.6. 

(c) Let Tc = May 1, 2008, be the time at which British Airways calculates the mark to 
market value of the CDS contract. We keep the CIR parameters of British Airways 
and Royal Dutch Shell at the same values calibrated in (a). We vary the volatility of 
the CIR process associated to Lehman Brothers, while keeping the other parameters 
fixed. We take the market CDS quotes of Lehman Brothers, British Airways and Royal 
Dutch Shell and recompute the shift process as 



^,(t;A) = log 



in > t) 



market 



in > t) 



market 



for any Tc < t < Td, where T^ = May 1, 2013. Here the market survival probabilities 
Qin > t)market, are Stripped from the CDS quotes of the three names at the date. 
May, 1, 2008, using a piece- wise linear hazard rate function. We compute the value 
CDSf^(rc, S'l, LGDi) of a risk adjusted CDS contract starting at and maturing at 
Td, where the five year running spread premium as well as the loss given defaults of 
the three parties are the same as in Tq. We have 

CDSfrf(Te, Si, LGDo,i,2) = CDS,,,(T„ Si, LGDi) - BR-CVA-CDS,,d(r„ Si, LGDo,i,2) 

(6.3) 

where CDSc,d(T'c, Si, LGDi) is the value of the equivalent CDS contract which does not 
account for counterparty risk and BR-CVA-CDSc,d(Tc, Si, LGDi) is the adjustment for 
the period [Tc, calculated at time Tc. 

(d) We calculate the mark-to-market value of the CDS contract as follows: 

MTMa^Si, LGDo,i,2) = CDSj^diTc, Si, LGDo,i,2) - ^^^a,bi^a,^^,^GDo,i,2) ^^^^^ 



Table 12 reports the MTM value of the CDS contract between British Airways and 
Lehman Brothers on default of Royal Dutch Shell under a number of correlation scenar- 
ios. The CDS contract agreed on January 5, 2006 is marked to market on May, 1, 2008 by 
British Airways using the four-step procedure described above. We check the effect of the 
increasing riskiness of Lehman Brothers by varying the volatility of the CIR process asso- 



ciated to Lehman. Table 12 reveals sensitivity of the BR-CVA to both default correlation 



and credit risk volatility of Lehman. We notice the following behavior. If British Airways 
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(roi,ro2,ri2) 


Vol. parameter ui 
CDS Impled vol 


0.01 
1.5% 


0.10 

15% 


0.20 

28% 


0.30 

37% 


0.40 

42% 


0.50 

42% 


(-0.3, -0.3, 0.6) 


(LEH Pay, BAB Rcc) 
(BAB Pay, LEH Rec) 


39.1(2.1) 
-84.2(0.0) 


44.7(2.0) 
-83.8(0.1) 


51.1(1.9) 
-83.5(0.1) 


58.4(1.4) 
-83.8(0.1) 


60.3(1.7) 
-83.8(0.2) 


63.8(1.1) 
-83.8(0.2) 


(-0.3, -0.3, 0.8) 


(LEH Pay, BAB Rcc) 
(BAB Pay, LEH Rcc) 


13.6(3.6) 
-84.2(0.0) 


22.6(3.2) 
-83.9(0.1) 


35.2(2.6) 
-83.6(0.1) 


43.7(2.0) 
-83.9(0.1) 


45.3(2.4) 
-83.9(0.2) 


52.0(1.4) 
-83.8(0.2) 


(0.6, -0.3, -0.2) 


(LEH Pay, BAB Rcc) 
(BAB Pay, LEH Rcc) 


83.1(0.0) 
-55.6(1.8) 


81.9(0.2) 
-58.7(1.7) 


81.6(0.3) 
-66.1(1.4) 


82.4(0.3) 
-71.3(1.1) 


82.6(0.3) 
-73.2(1.0) 


82.8(0.4) 
-74.1(0.9) 


(0.8, -0.3, -0.3) 


(LEH Pay, BAB Rec) 
(BAB Pay, LEH Rec) 


83.9(0.0) 
-36.4(3.3) 


82.9(0.1) 
-41.9(3.0) 


82.3(0.3) 
-55.9(2.2) 


82.9(0.2) 
-63.4(1.6) 


82.9(0.3) 
-65.8(1.5) 


83.0(0.3) 
-66.4(1.5) 


(0, 0, 0.5) 


(LEH Pay, BAB Rcc) 
(BAB Pay, LEH Rcc) 


50.6(1.5) 
-80.9(0.2) 


54.3(1.5) 
-80.5(0.3) 


59.2(1.5) 
-80.9(0.4) 


64.4(1.1) 
-82.3(0.3) 


65.5(1.3) 
-82.6(0.3) 


68.8(0.8) 
-82.8(0.3) 


(0, 0, 0.8) 


(LEH Pay, BAB Rec) 
(BAB Pay, LEH Rec) 


12.3(3.5) 
-80.9(0.2) 


21.0(3.0) 
-81.5(0.2) 


34.9(2.5) 
-81.9(0.3) 


41.3(2.1) 
-81.9(0.4) 


44.6(1.9) 
-82.1(0.4) 


50.6(1.4) 
-82.7(0.3) 


(0, 0, 0) 


(LEH Pay, BAB Rec) 
(BAB Pay, LEH Rec) 


78.1(0.2) 
-81.6(0.2) 


77.9(0.3) 
-81.9(0.2) 


79.5(0.5) 
-82.3(0.3) 


79.5(0.5) 
-82.2(0.4) 


80.1(0.6) 
-82.7(0.3) 


82.1(0.4) 
-83.2(0.3) 


(0, 0.7, 0) 


(LEH Pay, BAB Rec) 
(BAB Pay, LEH Rec) 


77.3(0.3) 
-81.2(0.2) 


77.3(0.4) 
-81.8(0.2) 


78.5(0.5) 
-81.9(0.3) 


79.2(0.5) 
-80.8(1.3) 


79.7(0.6) 
-82.4(0.3) 


81.5(0.4) 
-82.6(0.3) 


(0.3, 0.2, 0.6) 


(LEH Pay, BAB Rcc) 
(BAB Pay, LEH Rec) 


54.1(1.4) 
-81.3(0.2) 


56.7(1.3) 
-81.7(0.2) 


62.5(1.1) 
-81.4(0.4) 


63.6(1.1) 
-81.3(0.5) 


66.4(0.9) 
-81.6(0.4) 


69.7(0.6) 
-82.1(0.4) 


(0.3, 0.3, 0.8) 


(LEH Pay, BAB Rec) 
(BAB Pay, LEH Rcc) 


22.8(4.2) 
-83.0(0.2) 


28.8(3.5) 
-83.2(0.2) 


38.6(2.9) 
-82.8(0.3) 


42.6(2.9) 
-82.4(0.4) 


45.9(2.5) 
-82.5(0.4) 


52.0(2.2) 
-82.9(0.4) 


(0.5, 0.5, 0.5) 


(LEH Pay, BAB Rcc) 
(BAB Pay, LEH Rec) 


62.8(0.8) 
-67.4(1.1) 


64.5(0.8) 
-70.4(0.9) 


67.7(0.8) 
-72.9(0.9) 


68.5(0.9) 
-74.4(0.9) 


71.3(0.7) 
-75.8(0.8) 


73.2(0.6) 
-76.7(0.7) 


(0.7, 0, 0) 


(LEH Pay, BAB Rec) 
(BAB Pay, LEH Rcc) 


77.4(0.2) 
-47.3(2.2) 


77.3(0.3) 
-55.0(1.9) 


78.9(0.5) 
-61.6(1.6) 


79.1(0.5) 
-65.0(1.5) 


79.9(0.5) 
-67.5(1.3) 


81.4(0.4) 
-69.6(1.1) 



Table 12: Value of the CDS contract between British Airways and Lehman Brothers on 
default of Royal Dutch Shell agreed on January 5, 2006 and marked to market by Lehman 
Brothers on May 1, 2008. The pairs (LEH Pay, BAB Rec) and (BAB Pay, LEH Rec) denote 
respectively the mark-to-market value when Lehman Brothers is the CDS receiver and CDS 
payer. The mark-to-market value of the CDS contract without risk adjustment when Lehman 
Brothers is respectively payer (receiver) is 84. 2 (-84. 2) bps, due to the widening of the CDS 
spread curve of Royal Dutch Shell. 



is negatively correlated or uncorrelated with Royal Dutch Shell, see triples (0.6, —0.3, —0.2), 
(0.8,-0.3,-0.3), and (0.7,0,0), then the mark-to-market value of the risk adjusted CDS 
contract appears to be the largest for Lehman (this is true both if Lehman is the CDS payer 
and the CDS receiver). The adjustments are quite sensitive to the credit spreads volatility 
of Royal Dutch Shell. In particular, it appears from Table 12 that increases in credit spreads 
volatility of Shell increase the mark-to-market valuation of the CDS contract when Lehman 
is the CDS payer and decrease the contract valuation when Lehman is the CDS receiver. 
This is the case because larger credit spreads volatility increases the number of scenarios 
where the counterparty British Airways precedes Shell in defaulting. When Lehman is the 
CDS payer, this translates in larger CDS contract valuation for Lehman, as the negative 
adjustment done at the counterparty default time if the option on the residual net present 
value is in the money (and this is the the default intensity of Shell in 2008 has 

increased) does not take place. Conversely, when Lehman is the CDS receiver, this implies 
smaller CDS contract valuations for Lehman due to a symmetric reasoning. 

We next invert the role of Lehman Brothers and Royal Dutch Shell, hence Royal Dutch 
Shell becomes the investor and Lehman Brothers the reference entity. 

We consider the following two CDS contracts. In the first contract Royal Dutch Shell 
buys 5-year protection on Lehman from British Airways on January 5, 2006. In the second 
contract. Royal Dutch Shell sells 5- year protection on Lehman to British Airways on January 
5, 2006. As in the earlier case, British Airways computes the mark-to-market value of the 
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two CDS contracts on May 1, 2008. 



The results are reported in Table 13 



(roi,ro2,ri2) 


Vol. parameter ui 
CDS Impled vol 


0.01 
1.5% 


0.10 

15% 


0.20 

28% 


0.30 

37% 


0.40 

42% 


0.40 

42% 


(-0.3, -0.3, 0.6) 


(BAB Pay, RDSPLC Rec) 
(RDSPLC Pay, BAB Rec) 


513.0(1.9) 
-520.0(0.2) 


512.4(1.9) 
-520.0(0.2) 


512.8(1.9) 
-520.0(0.2) 


512.7(1.9) 
-520.0(0.2) 


513.4(1.9) 
-520.0(0.2) 


514.0(1.8) 
-520.1(0.2) 


(-0.3, -0.3, 0.8) 


(BAB Pay, RDSPLC Rec) 
(RDSPLC Pay, BAB Rec) 


511.1(2.3) 
-520.0(0.2) 


511.1(2.4) 
-520.0(0.2) 


511.4(2.3) 
-520.0(0.2) 


511.3(2.3) 
-520.0(0.2) 


511.9(2.3) 
-520.0(0.2) 


513.0(2.2) 
-520.1(0.2) 


(0.6, -0.3, -0.2) 


(BAB Pay, RDSPLC Rec) 
(RDSPLC Pay, BAB Rec) 


525.9(0.1) 
-442.2(3.6) 


525.9(0.1) 
-442.7(3.5) 


525.9(0.1) 
-442.6(3.5) 


525.9(0.1) 
-442.6(3.5) 


525.9(0.1) 
-442.3(3.5) 


525.8(0.1) 
-441.2(3.6) 


(0.8, -0.3, -0.3) 


(BAB Pay, RDSPLC Rec) 
(RDSPLC Pay, BAB Rec) 


526.5(0.0) 
-405.1(5.2) 


526.5(0.0) 
-405.3(5.2) 


526.5(0.0) 
-403.3(5.3) 


526.5(0.0) 
-405.4(5.2) 


526.5(0.0) 
-407.6(5.1) 


526.5(0.0) 
-403.9(5.3) 


(0, 0, 0.5) 


(BAB Pay, RDSPLC Rec) 
(RDSPLC Pay, BAB Rec) 


516.0(1.4) 
-503.7(0.8) 


515.7(1.4) 
-503.7(0.8) 


515.5(1.5) 
-503.7(0.8) 


516.3(1.3) 
-503.6(0.8) 


515.6(1.5) 
-503.8(0.8) 


516.7(1.3) 
-503.9(0.8) 


(0, 0, 0.8) 


(BAB Pay, RDSPLC Rec) 
(RDSPLC Pay, BAB Rec) 


511.6(2.4) 
-503.7(0.8) 


511.5(2.4) 
-503.7(0.8) 


512.1(2.2) 
-503.7(0.8) 


512.1(2.2) 
-507.5(4.0) 


505.7(7.1) 
-503.7(0.8) 


512.6(2.2) 
-503.9(0.8) 


(0, 0, 0) 


(BAB Pay, RDSPLC Rec) 
(RDSPLC Pay, BAB Rec) 


524.2(0.3) 
-504.2(0.8) 


524.2(0.3) 
-504.1(0.8) 


524.1(0.3) 
-504.2(0.8) 


524.2(0.3) 
-504.0(0.8) 


524.2(0.3) 
-504.2(0.8) 


524.2(0.3) 
-504.3(0.8) 


(0, 0.7, 0) 


(BAB Pay, RDSPLC Rec) 
(RDSPLC Pay, BAB Rec) 


524.8(0.3) 
-504.3(0.8) 


524.7(0.4) 
-504.3(0.8) 


524.8(0.3) 
-504.4(0.8) 


524.7(0.3) 
-504.2(0.8) 


524.7(0.3) 
-504.4(0.8) 


524.8(0.3) 
-504.5(0.8) 


(0.3, 0.2, 0.6) 


(BAB Pay, RDSPLC Rec) 
(RDSPLC Pay, BAB Rec) 


516.6(1.5) 
-484.3(1.7) 


517.0(1.5) 
-484.4(1.7) 


516.6(1.5) 
-484.3(1.7) 


517.3(1.4) 
-484.4(1.6) 


517.1(1.5) 
-484.5(1.6) 


517.2(1.4) 
-484.3(1.7) 


(0.3, 0.3, 0.8) 


(BAB Pay, RDSPLC Rec) 
(RDSPLC Pay, BAB Rec) 


507.4(5.5) 
-487.0(6.0) 


505.6(5.7) 
-484.5(1.7) 


508.9(4.5) 
-490.6(4.9) 


508.1(3.5) 
-492.7(9.4) 


502.6(7.7) 
-487.1(2.7) 


497.5(14.2) 
-488.3(3.3) 


(0.5, 0.5, 0.5) 


(BAB Pay, RDSPLC Rec) 
(RDSPLC Pay, BAB Rec) 


519.6(1.1) 
-460.2(2.8) 


519.6(1.1) 
-460.2(2.8) 


519.1(1.2) 
-459.6(2.8) 


519.7(1.1) 
-460.1(2.8) 


519.7(1.1) 
-459.4(2.8) 


519.4(1.1) 
-458.1(2.8) 


(0.7, 0, 0) 


(BAB Pay, RDSPLC Rec) 
(RDSPLC Pay, BAB Rec) 


523.8(0.3) 
-426.3(4.3) 


523.8(0.3) 
-426.0(4.3) 


523.7(0.3) 
-426.5(4.3) 


523.8(0.3) 
-427.5(4.2) 


523.8(0.3) 
-428.8(4.2) 


523.9(0.3) 
-424.4(4.4) 



Table 13: Value of the CDS contract between British Airways and Royal Dutch Shell on 
default of Lehman Brothers agreed on January 5, 2006 and marked to market by British 
Airways on May 1, 2008. The pairs (RDSPLC Pay, BAB Rec) and (BAB Pay, RDSPLC 
Rec) denote respectively the mark-to-market value when British Airways is the CDS receiver 
and CDS payer. The mark-to-market value of the CDS contract without risk adjustment 
when British Airways is respectively payer (receiver) is 529(-529) bps, due to the widening 
of the CDS spread curve of Lehman Brothers. 

The risk-adjusted mark-to-market value of the CDS contract is sensitive to correlation 
and it follows a pattern similar to the one discussed earlier for the other CDS contract. How- 
ever, differently from the previous case, there is less sensitivity to credit spreads volatility. As 
the default intensity of Lehman on May 1 2008 is already high, we have that increases in its 
credit spreads volatility does not vary significantly the number of scenarios where Lehman is 
the first to default. Therefore, the risk-adjusted mark-to-market value of the CDS contract 
does not depend much on credit spreads volatility in this case. 



7 Conclusions 

We have provided a general framework for calculating the bilateral counterparty credit val- 
uation adjustment (BR-CVA) for payoffs exchanged between two parties, an investor and 
her counterparty. We have then specialized our analysis to the case where also the under- 
lying portfolio is sensitive to a third credit event, and in particular to the case where the 
underlying portfolio is a credit default swap on a third entity. We have then developed a 
Monte-Carlo numerical scheme to evaluate the formula and thus compute the BR-CVA in 
the specific case of credit default swap contracts. We have provided a case study in Section 
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|5} and experimented with different levels of credit risk and credit risk volatilities of the three 
names as well as with different scenarios of default correlation. The results obtained confirm 
that the adjustment is sensitive to both default correlation and credit spreads volatility, 
having richly structured patterns that cannot be captured by rough multipliers. This points 
out that attempting to adapting the capital adequacy methodology (Basel II) to evaluating 
wrong way risk by means of rough multipliers is not feasible. Our analysis confirms that also 
in the bilateral-symmetric case, wrong way risk - namely the supplementary risk that one 
undergoes when the correlation assumes the worst possible value - has a structured pattern 
that cannot be captured by simple multiphers applied to the zero correlation case. 
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A Proof of the general counterparty risk pricing for- 
mula 

We next prove the proposition. 
Proof. We have that 

n(t,T) = C ASHPLOWs(t, T) 

= lyiU-B Cashflows (t, T) + l(7(j£)CASHFLOws(t, T) + li5(ji?CASHFLOws(t, T) (A.l) 

since the events in Eq. (2.4) form a complete set. From the linearity of the expectation, we 
can rewrite the right hand side of Eq. ( 2.6[ ) as 



Et{U{t, T) + LGDo ■ UuB -Dit, tq) ■ [-NPV(ro)]+ - LGD2 ■ lcuD-D{t, T2) ■ [NPV(r2)] + } (A.2) 
We can then rewrite the formula in Eq. ( |A.2 ) using Eq. A.l as 



E^l^UB Cashflows (t, T) + (1 - Rbc,o)1 AuBD{t,To)[-NPV{ToT 

+ lc7UDCASHFLOWs(t,T) + (ReC,2 - 1) IcuD^^ (^ , ^2) [A^PV^ (ts)] + 
+ 1 i5(ji?CASHFLOWs(t, T)] 

E4lAuBCASHFLOWs(t,T) + (1 - RBC,o)lAuBD{t,To)[-NPV{ToT] 

+ E4lcuDCASHFLOWs(t,T) + (ReC,2 - l)lcuD D{t , T2)[N PV {t^)]^] 

[l£;UJ7'CASHFLOWs(t, T)] (A. 3) 



We next develop each of the three expectations in the equality of Eq. (A. 3). 
The expression inside the first expectation can be rewritten as 

lAuBCASHFLows(t,T) + (1 - RBc,o)lAuBD{t,To)[-NPV{ro)] + 
= lAuBCASHFLows(t, T) + lAuBD{t, To) [~NPV{ro)]+ - RbcoIaubDH, Tq) [-NPV{to)] + 

(A.4) 
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Conditional on the information at tq, the expectation of the expression in Eq. (A. 4) is equal 
to 

Er, [l^uBCASHFLOWs(t,T) + lAuBD{t,To) {-NPV{To)f - Rbc, ol AUBD{t , To)[- N PV {to)]+] 

= [lyl(jij[CASHFLOWs(t, To) + Z)(t, ro)CASHFLOWs(To, T) +Z)(t,ro) (— [CaSHFLOWs(To, T)])^ 

-R^c.oD{t,To)[-NPV{ro)]^]] 

= l^(jB[CASHFLOWs(t, Tq) + D(t,TQ)'Kj-Q [CaSHFLOWs(To, T)] + D(t,T(j) (—]£,-(, [CaSHFLOWs(To, T)] ) 

-RFc,oD(t,ro)[-iVPV(ro)]+] 
= lAuB[CASHFLows(t,ro) + D{t, To) (E^o [CASHFLows(ro,T)])+ - REc,oD{t, to)[-NPV{to)]^] 
= l^uB[CASHFLows(t,ro) + Dit,To) (iVPF(ro))+ - REc,oZ}(t, ro)[-iVPy(ro)]+] 



where the first equality in Eq. (A. 5) follows because 

lAUi?CASHFLOWs(t, T) = l^yij [CASHFLOWs(t, Tq) + -D (t , Tq) CaSHFLOWS (Tq , T)] 



(A.5) 



being the default time Tq always smaller than T under the event AU B. Conditioning the 



obtained result on the information available at t, and using the fact that Et[E. 
due to t < To, we obtain that the first term in Eq. (A. 3) is given by 



TO [ 



[UuB [CASHFLows(t,ro) + D{t,To)iNPViTo))+ - REc,oD(t,ro)(-Arpy(ro))+]] (A.6) 
which coincides with the expectation of the third term in Eq. ( |2.5 ). 



We next repeat a similar argument for the second expectation in Eq. (A.3). We have 



lcuDCASHFLows(t,T) + (Rec.2 - 1 ) l^uD^ , ^2) [ A^P V (xa)] + 

= lcuDCASHFLOWs(t,T) - lcuDD{t,T2)[NPV{T2)]+ + R^C,2lcuDD{t,T2)[NPV{T2)] + 

(A.7) 

Conditional on the information available at time T2, we have 

E., [lcuDCASHFLOWs(t,T) - lcuDD{t,T2) {NPV{T2))-^RBC,2lcuDD{t,T2)[NPV{T2)]+] 
= Et-2 [lc'ur>[CASHFLOws(t, T2) + D(t, Cashflows , T) — -D(t, T2) (Et-2 [Cashflows (t2, T^)])"*" 

+R^C.2D{t,T2)[NPV{T2)]^]] 
= lcu/)[CASHFLOws(t, T2) + D(t, r2)ET-2 [CASHFLOws(r2, T)] — D(t, T2) (E^-j [Cashflows (r2, T)])"^ 

+R^c,2D{t,T2)[NPV{r2)]+] 

= ICUD [CASHFLOWs(t, T2) - D{t, T2) (E,2 [-CASHFLOWs(r2, r)])+ + REC,2D(t, r2 ) [A^P V (r2 )] +] 

= IcuD [CASHFLows(t, rs) - D{t,r2)i-NPV{T2))+ + REc,2Z}(t, r2)[ArPV(r2)]+] (A.S; 
where the first equality follows because 

lcuDCASHFLOws(t, T) = [Cashflows (t, T2) + Z)(i(:, r2)CASHFLOws(T2, T)] (^-9) 

being the default time T2 always smaller than T under the event C U D. Conditioning the 
obtained result on the information available at t < r2, we obtain that the second term in 
Eq. (A.S) is given by 

Et [icuD [CASHFLows(t, r2) + P)(t,r2)REc.2(Arp\/(r2))+ - P)(t, r2)(-Arp\/(r2))+]] (A.IO) 
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which coincides exactly with the expectation of the second term in Eq. (2.5). 



The third expectation in Eq. (A. 3) coincides with the first term in Eq. (2.5), therefore 
their expectations ought to be the same. Since we have proven that the expectation of each 
term in Eq. (2.5) equals the expectation of the corresponding term in Eq. (A. 3), the desired 
result is obtained. □ 



B Proof of the survival probability formula 

Proof. We have 

lcuDln>r.Q(ri >t|g.J 

= lr2<Tlr2<ro iU<T2<Ti + lr2<t lTi>r2lE [Q(Ai (t) < 6 I Sr2 , ^l) I Sr2] ) 
= lr2<Tlr2<rolA + ^[1 Q(Ai(t) <6|Sr2,6)l 

T2<t'^Tl>T2^[^Al(t)-Al{T2) 

(6-Ai(T-2))|S.2,{n>r2},{ro>r2}]) 

(B.l) 

The last step follows from the fact that the Ai{t) < is the same as Ai(t) — Ai(r2) < 
^1 — Ai(r2) and the right hand side — Ai(r2) becomes known once we condition on and 
St2- Here, by -FAi(t)-Ai(T2)) we indicate the cumulative distribution function of the integrated 
(shifted) CIR process Ai(t)—Ai(r2). Let us denote ?7ij = 1— exp(— Aj(r,)), where < i,j < 3 
denote the three names under consideration. Since = — log(l — Ui) and ti = A]'^(,^i), we 
can rewrite the inner term in Eq. ( B.l[ ) as 



1a + l.,<a.,>.,E[FA,(t)_Ai(r2)(-log(l - Ui) - Ai(r2))|g.„ {6 > Mr^)}, {^o > Ao(r2)}] 

= 1^ + 1^2<tl^,>^2E[FA,(t)_A,(^2)(-log(l - Ui) - Ai{t2))\St2,Ui > f7l,2,f/o > Uo,2] 

(B.2) 



Then we can rewrite Eq. (B.2) as 
"1 



1^ + 1 



T2<t-^T1>T2 



^AiW-Ai(r2)(- log(l - ui) - Ai{r2))dQ{Ui < ni|g,,, t/i > t/1,2, f/o > Uo,2) 

;-log(l - Ml) - Ai(r2))rfQ([/i < Ui\5r2,Ui > f7i,2, Uo > 1/0,2) 

(B.3) 



T2<t-I-T1>T2 l_ -rAl(t)-Al(T2) 
Ul,2 



The conditional distribution may be computed as follows. Denote 



Ci\o,2{ui; U2) 



< ui\%,,Ui > Uio,Uo > f/0,2) 



(B.4) 



Eq. (B.4) may be rewritten as 



Cl\o,2{ui-U2) = Q{Ui < m\U2, Ui>U 1,2, Uo>U 0,2) 

QiU^ < ni, U^> U,,2\U2, Uo_> ^0,2) 

QiU,>U,,2\U2,Uo>Uo,2) 
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Notice that U2 is explicitly known once we condition on and is given by 1 — exp(A2('r2)). 



Let us next express Eq. (B.5) in terms of the copula function. We start with the numerator, 



which may be expressed as 

h < Ml, Ui > Ui^2\U2, Uo > 1/0,2) 



Q(f/i <ui,Ui>Ui,2,Uo>Uo,2 


U2) 


Q(f/0>f70,2|f/2) 




Q{U,<ui,Uo>UoMU2)-QiL 


h<Ui,2,Uo>Uo,2\U2) 


Q(f/0 > Uo,2 


U2) 



(B.6) 



We have that 



1 < Ui,Uo > f/o,2|f/2) 



h < ui\U2 

dCi^2iui,U2 



1 < Mi,f/o < f/o,2|f/2) 

dC{U 0,2, ui,U2) 



U2=U2 



(B.7) 



U2 = U2 



where Ci^2 denotes the bivariate copula connecting the default times of names "1" and "2" , 
while C denotes the trivariate copula. Similarly, we have that 



1 < Ui^2,Uo > Uo^2\U2) 



dCi^2iUi^2,U2) 



dC{Uo,2,U^,2,U2) 



U2=U2 



(B.8) 



U2=U2 



The denominator in Eq. (B.5) may be computed as 



1 > U^2\U2,Uo > Uqo) 



where 

h>Ul,2,Uo>Uoo\U2) 



> Uo2\U2 



^ _ dCo^2{,Uo,2,U2) 
dU2 



1 > Uio.Uq > Uoo\U2 



> Uo2\U2 



> Uo,2,Ui < t/i,2|f/2 
_ 9Cl,2(^7l,2,M2 

U2 = U2 



(B.9) 



+ 



U2 = U2 



dC{Uo,2,Ul,2,U2) 

duo 



U2=U2 

(B.IO) 



All together, we have that 

9Cl,2(ltl,M2) 

Ci|o,2(mi; f/2) = 



du2 



dC{Uo^2,ui,U2) 



U2 = U2 



du2 



dCl^2{Ul,2,U2) 



U2 = U2 



du2 



+ 



aC(C/0.2,C/l.2,»2) 



tt2 = f^2 



du2 



U2=U2 



i9C(),2(C/o,2,M2) 

du2 



U2=U2 



dC 1 ^2 (U 1,2,112) 
du2 



+ 



U2=U2 



dCiUo,2,Ul,2,U2) 

du2 



U2 = U2 



(B.ll) 



□ 



The case when the investor is the first to default is symmetric, leading to the conditional 
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distribution given by 
Ci\2,o{ui]Uo) := 



1 < Ui\U2, {n > To}, {t2 > To}) 

1 < ui\Uo, Ui > Uifl, Uo > U2,o) 



dCo^i{uo,ui) 



duo 



dC{uo,ui,U2.o) 



uo=Uo 



duo 



dCo^i{uo,Ui,o) 



uo=Uo 



duo 



+ 



dC{uo,Ul,o,U2,o) 



uo=Uo 



duo 



uo=Uo 



1 - 



duo 



dCo^i{uo,Uifi) 



uo=Uo 



duo 



+ 



dC{uo,Ui,o,U2fi) 



uo=Uo 



duo 



uo=Uo 



(B.12) 



and the survival probability conditioned on the information known at the investor default 
time is given by 



lAuBlTi>ToQ(n >t\S 

'^to<t'^to<T2 



Ifi + lro<tlri>ro /_ ^Ai(i)-Ai(ro) log(l " Ml) " Al (Tq) )rfCi|2,0 («! ; f^o) 
-'C/1,0 / 



where 



B = {t<To< Ti} 



(B.13) 
(B.14) 
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